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We address the question of the multiplicativity of the maximal p-norm output purities of bosonic Gaussian 
channels under Gaussian inputs. We focus on general Gaussian channels resulting from the reduction of unitary 
dynamics in larger Hilbert spaces. It is shown that the maximal output purity of tensor products of single-mode 
channels under Gaussian inputs is multiplicative for any p £ (1, oo) for products of arbitrary identical channels 
as well as for a large class of products of different channels. In the case of p = 2 multiplicativity is shown to be 
true for arbitrary products of generic channels acting on any number of modes. 
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I. INTRODUCTION 

Additivity and multiplicativity questions play a central role 
in the field of quantum information theory: does it help to 
make joint use of a quantum channel when transmitting quan- 
tum or classical information in form of entangled inputs, or 
is one better off by merely invoking the channel many times 
with uncorrelated inputs? Or, what is the entanglement cost, 
the rate at which maximally entangled pairs need to be in- 
vested in the asymptotic preparation of a mixed bi-partite state 
using local operations and classical communication only? If 
the so-called entanglement of formation turned out to be addi- 
tive, the evaluation of this quantity, which amounts to a much 
simpler optimization problem, would be sufficient to provide 
the complete answer to this question. Most instances of such 
additivity problems in quantum information share the com- 
mon feature of being notoriously difficult to solve. Recently, 
yet, a picture emerged that made clear that several of these 
problems share more than a formal similarity Ql |2> S 01 ■ It 
has actually been shown that at least four instances of such 
additivity problems are logically equivalent, being either all 
wrong or all true. Besides the fundamental insight that this 
equivalence provides, such an observation is practically help- 
ful, since it links isolated additivity results to other instances 
of such problems. 

This paper is concerned with a specific multiplicativity 
question for quantum channels of bosonic systems: it deals 
with the maximal output purity of Gaussian channels under 
Gaussian inputs. This quantity specifies how well the purity 
of an input state, measured in terms of p-norm purities (or 
equivalently Renyi entropies), can be preserved under the ap- 
plication of the (generally decohering) channel and provides a 
way to characterize the decoherence rate of the channel. If this 
output purity turns out to be multiplicative for a tensor product 
of channels, then input entanglement cannot help to better pre- 
serve the coherence of the output states jp|]. The multiplica- 
tivity of the maximal output purity for p ^ 1+ corresponds 
to the additivity of the minimal von Neumann entropy. Fur- 
thermore, if general inputs are allowed for, the multiplicativity 



for such a limiting instance is strictly related to the additivity 
of the (appropriately constrained) Holevo capacity and of the 
entanglement of formation (EoF) ll|, 0,131 ■ 

In turn, the quantum information properties of bosonic 
Gaussian channels Q 01 have attracted strong theoretical at- 
tention in recent years. This class of quantum channels is 
practically very important: indeed, the transmission of light 
through a fiber is described to a very good approximation by a 
Gaussian bosonic channel. The unavoidable coupling to exter- 
nal field modes yield losses, whereas excess noise can be in- 
corporated as random classical Gaussian noise, reflecting ran- 
dom displacements in phase space. The estimation of various 
information capacities, both quantum and classical, has been 
thoroughly addressed for these Gaussian channels H. As for 
the multiplicativity of the maximal output purity, the quest has 
been challenged in a series of recent works, and strong argu- 
ments have been provided to support it, addressing a subset 
of channels investigated in the present paper j^l- In fact, for 
a specific channel model describing a beam splitter interac- 
tion of a bosonic mode with a thermal noise source and for 
integer p > 2 multiplicativity of the maximal p-norm output 
purity was recently proven f^. However, a definitive proof of 
the multiplicativity conjecture for tensor products of general 
Gaussian channels and non-integer p is still missing. 

In the present paper we aim to make a step towards a theory 
of channel capacities of Gaussian channels under Gaussian in- 
puts, dealing with more general instances of channels. Such 
a setting, besides being interesting in its own right, yields ob- 
vious bounds for the unconstrained maximal output purity for 
Gaussian bosonic channels, which in turn may be conjectured 
to be tight as it is the case for particular single-mode channels 
and integer p > 2 J31. Moreover, the present paper is meant to 
be a further step towards a clear picture of a general quantum 
information theory of Gaussian states, linking channel capac- 
ities with entanglement properties. This picture could provide 
a powerful laboratory, when a complete solution to the speci- 
fied additivity problems is lacking. 

The paper is structured as follows. In section Hll we intro- 
duce the notation, basic facts about Gaussian states and we 



define the class of Gaussian channels we will deal with. In 
section HII] the p-norms as measures of purity are presented 
and determined for Gaussian states. In section ITVl the Gaus- 
sian multiplicativity of the maximal output purity is defined, 
while section|y]contains all the analytical results about mul- 
tiplicativity. Finally, in section IVll we review our results and 
provide some comments and perspectives. 



II. GAUSSIAN STATES AND CHANNELS 

A. Gaussian states 

We consider quantum systems with n canonical degrees of 
freedom, i.e., a system consisting of n modes. The canonical 
coordinates corresponding to position and momentum will be 
denoted as i? = (ii,pi, ...,a;„,p„), where in terms of the 
usual creation and annihilation operators we have that x, 



(flj + a|)/y2 and 



Pi 



-i{ai — a|)/%/2. In terms of the 



Weyl operators 

the canonical commutation relations (CCR) can be written as 

Wl W^' = Wl W^> e*«^ '"«' , (2) 

where 



-® -1 



1 




(3) 



The latter matrix a is the symplectic matrix. States can be 
fully characterized by functions in phase space (R^", a). The 
characteristic function is defined as 



where the state p can in turn be expressed as 
1 



P = 



(2^)' 



rf'"?xp(-0w^«; 



(4) 



(5) 



The characteristic function is the ordinary Fourier transform 
of the Wigner function commonly employed in the phase 
space description of quantum optics llOtlllL 

Gaussian states are, by definition, the states with Gaussian 
characteristic function (and therefore Gaussian Wigner func- 
tion) 



X(0 = X(0)c 



-^ire+iD 



(6) 



Gaussian states are fully determined by first and second mo- 
ments of the quadrature operators, respectively embodied by 
the vector d = aD and by the real symmetric 2n x 2?! matrix 
r = a^^a, with 

lij — '^{^i^j + RjRi)p ^ {Ri)p{Rj)p J (7) 



where {0)p = tr [pO] for the operator O. First moments can 
be set to zero by local unitaries, so that they play no direct 
role in properties related to entanglement and mixedness of 
Gaussian states. To our aims a Gaussian state will be charac- 
terized by its covariance matrix 7. The covariance matrix 7 of 
a Gaussian state p (and, indeed, any covariance matrix related 
to a physical state), has to satisfy the uncertainty principle 



7 + icr > , 



(8) 



reflecting the positivity of p. Subsequently, Q will stand for 
the set of Gaussian states with vanishing first moments (for 
simplicity, the underlying number of canonical degrees of 
freedom will not be made explicit). Pure Gaussian states are 
those for which dot 7 = 1. These are the minimal uncertainty 
states, saturating Ineq. (|8}. The subset of pure Gaussian states 
with vanishing first moments will be denoted as Q. 

Any unitary U generated by polynomials of degree two in 
the canonical coordinates is, by virtue of the Stone-von Neu- 
mann theorem, the metaplectical representation of a real sym- 
plectic transformation S e Sp{2n, R). We recall that the real 
symplectic group Sp{2n, R) consist of those real 2n x 2n 
matrices S for which S^aS = a. Such symplectic opera- 
tions preserve the Gaussian character of the states and act by 
congruence on covariance matrices 



7 



S^jS 



On Weyl operators such an operation is reflected by 



(9) 



(10) 



We mention that ideal beam splitters and squeezers are de- 
scribed by symplectic transformations. The expression of the 
generators of Sp{2n, R) will be useful in the following and is 
detailed in App. |A| Moreover, we recall that a useful way to 
express a generic symplectic transformation S is provided by 
the Euler decomposition 1121 



S = O'ZO", 



(11) 



where O' ,0" G K{n) = S'p(2n, R) n SO{2n) are orthog- 
onal symplectic transformations, whose set forms the maxi- 
mal compact subgroup of Sp{2n, R). They are those oper- 
ations typically referred to as being passive, again, in opti- 
cal systems corresponding to beam splitters and phase shifts. 
The group of all Z — diag (zi, 1/zi, . . . , 2„, l/z„) with 
zi, ..., Zn € R"'"\{0} is the non-compact group of all such 
Z, reflecting local squeezings; this group will be denoted by 
Z{n) in the following. 

A frequently used tool will be the fact that any CM 7 can 
be brought to the Williamson normal form I13ll 



sis- - e 



k/ 

i/i 



(12) 



with i/l e [1, 00) and S e Sp{2n, R). The vector (i/|, ..., i/,\) 
is the vector of decreasingly ordered symplectic eigenvalues, 
which can be computed as the spectrum of the matrix \icr^\. 



The previous decomposition is nothing but the normal mode 
decomposition. Choosing the standard number basis {\n) : 
n G N} of the Hilbert space associated with each mode, the 
Gaussian state with vanishing first moments and the second 
moments as in the right hand side of Eq. ( I12> is given by 



P = 



^ fc=0 



-1 



lfc)(fcl. 



(13) 



Recalhng the Euler decomposition and Williamson's theo- 
rem, expressed by Eq. Jll> and Eq. (I12> one then finds that 
the CM 7 of an arbitrary pure Gaussian state reads 



7 = O^ZO , with O e K{n) , Z e Z{n) . 



(14) 



We finally mention that, as it is evident from the definition of 
the characteristic function, tensor products of Hilbert spaces 
correspond to direct sums of phase spaces. Therefore an un- 
corrected tensor product of Gaussian states with CMs 7^, 
i ~ 1, ...,n, has the CM 7 = (Bf^iJi- Likewise, for a 'lo- 
cal' tensor product of symplectic transformations Si one has 



B. Gaussian channels 



This class of channels includes the classical case of random 
displacements with a Gaussian weight 



P' 



j <f''^P{i)wlpW^, (20) 



where P(^) = P(0)e~25^ ? is a multivariate Gaussian 
with positive covariance matrix Y . In our notation, such a 
channel corresponds to X = 1, y > 0. The amplification 
and attentuation channels can be described by Eq. M%\ too, 
with 



X=el2, Y 



II 



e'\l. 



(21) 



with e < 1 (attenuation) or e > 1 (amplification). Note that 
this instance encompasses the case of white noise as well. 

Also the dissipation in Gaussian reservoirs after a time t is 
included in such a class of channels, with the choices 



X 



-rt/2 



1, y = (i 



-rt 



)lB. 



(22) 



where F is the coupling to the bath and 73 is the covariance 
matrix describing the reservoir 

The channel model on which the papers of Ref. |21 are 
focused on is characterized by X = cl, c < 1 and Y diagonal 
(mainly Y (x 1). 



In general, a Gaussian channel is a trace-preserving com- 
pletely positive map that maps Gaussian trace-class operators 
onto Gaussian trace-class operators. A Gaussian channel is 
defined by its action on the Weyl operators, according to 



Wi 



Wx(,e 



-y{i) 



(15) 



where X is a real 2n x 2n-matrix, and y is a quadratic form. 
We do not consider Unear terms in this quadratic form, which 
would merely correspond to a displacement, i.e., a change in 
first moments. We can hence write y{^) = ^^Y^/2. Com- 
plete positivity of the channel requires that 



Y 



iX^aX > 



(16) 



For single-mode channels, this requirement is equivalent to 



r > 0, det[y] > (det[X] - 1)^ 



(17) 



The second moments are transformed under the appUcation of 
such a channel according to 



7 



X^jX + Y 



(18) 



Any channel of the form of Eq. ( I18> corresponds to the re- 
duction of a symplectic (unitary) evolution acting on a larger 
Hilbert space and, vice versa, any evolution of this kind is de- 
scribed by Eq. ( I18> for some X and Y ilm . 

In the Schrodinger picture, we will denote such channels 
(characterized by the matrices X and Y) by ^x,y, acting as 



*-v.y(p) • 



(19) 



The additional noise term F of a general channel incorpo- 
rates both the noise that is due to the Heisenberg uncertainty 
in a dilation, and the additional classical noise. In the same 
manner as minimal uncertainty Gaussian states can be intro- 
duced, pure channels can be considered, satisfying 



Y = -{X^aX - a)Y-^{X^(jX - a), 



(23) 



where the inverse has to be understood as the Moore-Penrose 
inverse. 



III. MEASURES OF PURITY 



Generally, the degree of purity of a quantum state p can be 
characterized by its Schatten p-norm I14ll 

||p||p = (tr|pr)? = (trpP)?, pG(l,oo). (24) 

We mention that the case p = 2 is directly related to the quan- 
tity often referred to as linear entropy or purity in the closer 
sense, p = trp^ = WpW^. The p-norms are multiplicative 
on tensor product states and determine the family of Renyi 
entropies Sp H5i] . given by 



In tr pP 



(25) 



quantifying the degree of mixedness of the state p. It can be 
easily shown that 

lim Sp = — tr [p In p] = Sy (p) . (26) 



Thus the von-Neumann entropy Sy is determined by p-norms, 
as it is given by the first derivative of \\p\\p at p ^ 1+. The 
quantities Sp are additive on tensor product states. It is easily 
seen that Sp{p) E (0, oo), taking the value exactly on pure 
states. 

Because of the unitary invariance of the p-norms, the quan- 
tities tr pP of a n-mode Gaussian state p can be simply com- 
puted in terms of its symplectic eigenvalues. In fact, due to 
Eq. ( I12> . tr pP can be computed exploiting the diagonal state 
1/ of Eq. jl3> . One obtains 



tr pP = Y[ 



2P 



i^pn 



where 



1 fpiA) ^pW 



fp[x) = {x + iY~{x~iy 



(27) 



(28) 



and we have defined Fp{'-f) ~ n"=i fvi^i)^ ™ terms of the 
symplectic eigenvalues of the co variance matrix 7 of p. A first 
consequence of Eq. (I27t is that the purity /^j of a Gaussian state 
is fully determined by the symplectic invariant dct 7 alone: 



m(p) 



1 



n 



,i 



i=i ' 



v/det[7] 



(29) 



Every Renyi entropy, or p-norm respectively, yields an order 
within the set of density operators with respect to the purity of 
the states. Yet a stronger condition for one state being more 
ordered than another one is given by the majorization relation 
which gives rise to a half-ordering in state space. A density 
operator p is said to majorize p, i.e., p> p\f 



E^i^E^ 



(30) 



i=l 



for all r > 1, where AJ is the decreasingly ordered spectrum 
of p. Majorization is the strongest ordering relation in the 
sense that if p y p, then tr/(p) < tr/(p) holds for any con- 
cave function / and in particular for every Renyi entropy J14ll . 
It has recently been conjectured |9|| for a special class of Gaus- 
sian single-mode channels that the maximal output purity is 
not only achieved for a Gaussian input state but that the opti- 
mal Gaussian output even majorizes any other possible output 
state. 



Thus, on the level of second moments, the multiplicativity of 
the Gaussian maximal output p-purity under Gaussian inputs 
corresponds to the multiplicativity of the infimum of Fp 






p[(PX,,Yi 



(7)) , (33) 



where the infimum is taken over all covariance matrices. 

For finite dimensional systems and the usual definition of 
maximal output purity (allowing for input on the whole con- 
vex set of trace class operators), the convexity of the p-norms 
guarantees that the sup of Eq. ( I31> can be approached restrict- 
ing to pure states. The set Q of Gaussian states is not convex. 
However, every Gaussian state still has a convex decomposi- 
tion into pure Gaussian states such that it is again sufficient to 
consider pure input states only: 

Lemma 1. - For any Gaussian channel ^x,y ■ Q ^ Q and 
any p G (1, 00), one has 



sup 

peg 



\'Px,Y{p)\\p = snp\\^x.Yip)\\p. (34) 

peg 



Proof. Consider the Williamson standard form of the covari- 
ance matrix 7 = S^vS. By rewriting this as 



7 = 5^5 + S'^iv - 1)5 =:ip + V 



(35) 



one infers from Eq. ( I20> that the state corresponding to 7 can 
be generated by randomly displacing a pure state with covari- 
ance matrix 7p in phase space according to a classical Gaus- 
sian probability distribution with covariance V . Hence, the 
state has a convex decomposition into pure Gaussian states 
and the Lemma follows from the convexity of the p-norms.D 

Let us now consider a channel ^x,y resulting from the ten- 
sor product of the channels ^Xi.Yi ,i — I, ■■■,n, 



*x,Y =0*x.,F., 



(36) 



acting on Gaussian states associated with the tensor product 
Hilbert space. Since tensor products in Hilbert spaces corre- 
spond to direct sums in phase space, we have that ^x,y = 
^®Xi,®Yi- We will say that the Gaussian maximal output p- 
purity of the channel ^x.y is multiplicative if 



IV. MULTIPLICATIVITY OF THE MAXIMAL OUTPUT 
PURITIES 

We define now forp € (1, c)o) the Gaussian maximal output 
p-purity of a Gaussian channel ^x.y : 5 — > ^, as 



Cp(*x,y) =sup \\^x,y{p)\\p 
peg 



(31) 



where the sup is taken over the set of Gaussian states. In terms 
of covariance matrices and of the function Fn, one has 



2" 



U'^x,, 



p. 



MFp {(f>x,Y{l)) 
p 



(32) 



^pi^x.y) = Y[Cp{^x,.,Y). 



(37) 



j=i 



Let us remark that this is equivalent to stating that the max- 
imal output purity can be attained by means of uncorrected 
input states. More precisely, for tensor products of channels 
one has that, denoting by S the set of product Gaussian states 
with respect to each of the modes, the multiplicativity of the 
maximal output p-purity of the channel ^x,y is equivalent to 



^p(*x,y) = inf I 
pes 



*x,y(p)llp. 



(38) 



V. MULTIPLICATIVITY STATEMENTS 

From now on, we will mainly restrict to tensor products of 
single-mode channels, for which the matrices describing the 
channels are direct sums of 2 x 2 matrices Xi and Yi, i = 
1, ..., n: X = ®i=iXi and Y = ®f^iYi. Moreover, we will 
assume that the determinants of the Xi have equal sign. In this 
case the invariance of the p-norms and of the correlations of 
quantum states under local unitary operations can be exploited 
to simplify the problem, according to the following. 

Lemma 2. - Let the determinants of the matrices {Xi : 
1, ...,n} have equal signs. Then the Gaussian maximal out- 
put p-purity of the tensor product of single-mode channels 
^®Xi,mYi is multiplicative if and only if the Gaussian max- 
imal output p-purity of the channel <^ ^ y is multiplicative, 
with 



X = 01,, Y = ^%, (39) 

X, = V|det[X,]|l2, Y, = v/det[K,]l2. (40) 



Proof. Let us first reduce the case of negative determinants to 
that of dct[Xj] > 0. To this end we write Xi = azX^ , so 
that dct[X+] = - det[Xi]. Since 6I76' with 9 = 0"^^ a, is 
again an admissible covariance matrix (corresponding to the 
time reversed state) and in addition 9^ = 1, we have indeed 
that 



ini FpiX'-fX 

1 



Y) ^ inf Fp(X+'^7X+ + Y) . (41) 



Now let us see how the case det[Xi] > can be reduced 
to the standard form given in the Lemma. Due to the unitary 
invariance of the p-norm we can replace X^Yhy X ^ S'XS, 
Y = S'^YS, with S, S' being any symplectic transformations 
and 



[niFp{X'^-/X + y) = inf Fp{X'^-fX + Y) . 



(42) 



In particular we may choose S = 0'jLi SiOi, such that 
Oi e SO{2) and S*, e 5p(2,R) bring Y in standard form: 
SiYSf = Vdet Yi 1 = Yi. Furthermore, we choose S" to 
consist of blocks S[ = ZiO'i, Oi G SO{2), Z, e Z{2) such 
that the orthogonal matrices Oi, O'i diagonalize XiSi in 

X, = Z, [0[{X^Si)0^] (43) 

and the squeezing transformation Zi gives rise to equal diag- 
onal entries yielding Xi = y^dct[Xi]|l. According to the 
block structure of the involved transformations (correspond- 
ing to the direct sum of 'local' single mode operations) one 
has 

n n 

[| inf F.p{Xj-iX, + Y,) = \{ inf Fp(Xf 7X, + Y) . (44) 

i=l i=l 

Eqs. J42> and J44> straightforwardly imply that 

inf^ Fp{XT-iX + r) = nr=i inf 7 Fp{Xf-iX, + Yi){f and 



onlyifinf^Fp(X^7X + y) = nr=i inf7 ^p(^f 7^^ + f.)- 
This proves the claimed equivalence of the multiplicativity 
statements. D 

One remark is in order concerning the case det Xi = in 
the above Lemma. In fact, in this case multiplicativity is trivial 
and the maximal p-norm output purity does not at all depend 
on the Xi. To see this note that for two positive matrices A > 
S > we have i^i{A) > vl^B) |T3], which implies that 
inf^ Fp{X^X'^ + y) > Fp{Y)- This becomes, however, an 
equality in the case dct Xi = since we can always choose 
the input state to be a product of squeezed states such that in 
the limit of infinite squeezing X^X"^ -^ 0. 

A first relevant consequence of Lemma 2 follows. 

Proposition 1. - The Gaussian maximal output p-purity of a 
tensor product of n identical single-mode Gaussian channels 
^®iX,®iY is multiplicative for any p € (1,00). Moreover, the 
output corresponding to the optimal product input majorizes 
any other Gaussian output state of the channel. 

Proof. We recall that, because of Euler decomposition, the 
covariance matrix of any pure Gaussian state p can be written 
as 7 = O'^ZO, where O e K{n) = S'p(2n,]R) n SO{2n) 
is an orthogonal symplectic transformation and Z e Z{n) 
corresponds to a tensor product of local squeezings, Z — 
diag (zi, l/zi, . . . , z„, 1/ Zn). Clearly, if O = 1 then the state 
is uncorrected. For a tensor product of identical channels, 
Eq. (|^J holds globally, X = xl2n, Y = y\2n- Therefore, 
exploiting Lemma 2 and the invariance of Fp under symplec- 
tic transformations, one has 

ini FpiX-iX + Y)^ inf F^ix^O'^ ZO + yl) (45) 
= inf Fpl-x^Z + yl). (46) 

ZeZ(n) 

Due to the block structure of elements in Z{n) this proves the 
first part of the proposition. 

For the majorization part we exploit the fact that a compo- 
nentwise inequality for the symplectic eigenvalues vj < i)\ 
for all i impUes majorization on the level of density operators, 
i.e. p>- p ITtIi . The symplectic eigenvalue vj of the output co- 
variance matrix 7' = X'^^X + y is given by the square root 
of the ordinary eigenvalue AJ of the matrix a^' a^ ^' (when 
appropriately taking degeneracies into account). Continuing 
with the expression in Eq. ( I46> we have thus to consider the 
dependence of 



Ki^io'i) 



y^ + x^y\\{Z + Z-^) (47) 



on Z. However, choosing Z = 1 in Eq. (l47l minimizes all 
the eigenvalues simultaneously and thus proves the desired in- 
equalities between the optimal and any other Gaussian output 

state. D 



In the following we will investigate the multiplicativity is- 
sue in the case of tensor products of different Gaussian chan- 
nels. To proceed in this direction, we aim to turn our opti- 
mization problem over the non-convex set of Gaussian states 
into an analytical one. To do so two simple remarks, giving 



rise to alternative parametrizations of pure covariance matri- 
ces, will be exploited. Firstly, because of the Euler decompo- 
sition given by Eq. J14> . the set of the pure n-mode covariance 
matrices can be parametrized by means of the functions 



7 : (R+)" X R" 



defined as 



7(Z,z) = e-S-i'--^^£l(2)e-S"=i''^' 



(48) 



(49) 



where the Li, i = l,...,n?, are the generators of 
the compact subgroup K{n) (see App. 13 and D{z) ~ 
diag (zi, l/zi, . . . ,Zn, 1/zn) e Z{n), with Zi > for all i. 
Here, I = (Zi, ..., ^„2) is a vector of n^ real parameters while 
z = (zi, ..., z„) is a vector of n real strictly positive parame- 
ters. 

Otherwise, the set of pure covariance matrices admits the fol- 
lowing parametrization 



%{k) 



y' , kiK"^ y^ , kiKi 



(50) 



where 7 is an arbitrary pure covariance matrix, the Ki are the 
d = 2n^ + n generators of the symplectic group (detailed in 
App.lAl and k = (fci, . . . , ka) is areal vector of dimension d. 
We are now in a position to prove our main result. 

Proposition 2. - The Gaussian maximal output p-purity of a 
tensor product of single-mode Gaussian channels '^®Xi,®Yi 
with y^ > and identical Aet[Xi] for all i is multiplicative 

for any p G (1, 00). 

Proof. Because of Lemma 2, this multiplicativity issue is 
equivalent to the one for the 'simplified channel' $^ y, with 

X ^ a;l2„, X =_v/|det[X,]| and Y = ®,^Aq\}^,\\^ ac- 
cording to Eqs. ( 1391 140> . For ease of notation, and since the 
subsequent argumentation does not depend on the value of .t, 
we will state the proof for x = \. 

In a first step we aim to show that the infimum of Eq. ( I33> 
is indeed a minimum, that is, the infimum of Fp($j y(7)) 
is achieved for a defined input and not asymptotically ap- 
proached in the non compact set of pure covariance matrices. 
Therefore, we will analyse the asymptotic behaviour of the 
output purity of the channel ^^ y in the limiting case of infi- 
nite squeezing. For a given channel of this kind, let us define 
the function G^ y : (IR+)" x R"' ^ IR as 



G„ya,z) = (^pO$jy)(7(/,z)), p>l, 



(51) 



where the function ^(l^z) has been defined in Eq. ( I49t . 
To show that the function is indeed attained by a (possibly 
not unique) given covariance matrix, we address the asymp- 
totic behaviour of G y, showing that its infimum cannot be 
asymptotically approached. To see this, let us investigate the 
product of the symplectic eigenvalues 



nK^)' = det[7a,z) + y] 



(52) 



i=l 



Such a function is periodic in I, because these variables are 
related to the compact subgroup of Sp{2n, R) (consisting of 



rotations). Thus the domain of the I £ R" can be chosen 
compact. Therefore, only the cases Zi -^ and Zi ^ 00 have 
to be considered, for i = 1 , . . . , n. Let us take the ordered 

Ust of ?/|, ...,yJj, where Y = ©f^^ 



Yi andf, = yfla. Then, 



clearly, by the positivity of Y (implied by the positivity of Y), 
yl^ > 0. Now, it can be easily shown 1 14] that for two positive 
m X m matrices A and B with ordered lists of eigenvalues 



n, ■ 



i^i. and bi, ..., fe^i respectively, one has 

ra 

dct[A + B]>l[albl. 



(53) 



i=l 



Therefore, noticing that the values Zi , 1/zi constitute the spec- 
trum of 7(/, z), we get 

n ^ 

det[7(/, z) + Y]> [](z, + yi)(- + y^) . (54) 



From this it immediately follows that, for any i = 1, . 



lim dct[7(/,z) + y] 



lira dct[7(Z,z) 



Y] 



(55) 

Eq. (I55> shows that the product of the sympletic eigenvalues 
of 7(/, z) + Y diverges for Zi -^ and Zi -^ 00. Moreover, 
all symplectic eigenvalues of ^{l, z) + Y are clearly positive. 
This shows that the function G y diverges for Zi — > and 
Zi ^ 00 for alH = 1, ..., n. In turn, this means that the global 
infimum of the function Fp o ^-^ y-(7) over the set of pure 
covariance matrices is indeed a minimum. 

In what follows we show that such a minimum is achieved 
for 7=1, that is for a manifestely uncorrected input, thus 
completing the proof. Let us first notice that, since Fp is a 
function of the symplectic eigenvalues. 



FpiV^jV) = Fpi^) 



(56) 



for all^ e G,V e Sp{2n, R), and all channels $j y . 

Now, let us call M the infimum of the function Fp o $ ^ y on 
the whole space of covariance matrices. Suppose that 7 is one 
of the optimal matrices, granting the minimum of FpO (^^y. 
Let us denote by y a symplectic transformation bringing 7 + 



Y in Williamson form, so that V'^jV 
Exploiting Eq. ( I56t . we have 



V'^YV = ez/Ms 



M = Fp{^ + Y)^Fp(y^jV + V^YV). (57) 

For the sake of simplicity, let us define 7' = V'^^V and 
Y' ~ V^YV . We may write the matrix 7' in terms of 2 x 2 
submatrices as 



/ ai Pi 



PI2 



Pin \ 



Pn^ 



\Pl 



Pn-lr 



(58) 



For a given channel and covariance matrix 7, let us define 
the function Hp ^ Fp o ^i,Y'{%'{k)) : W^ — > R (with 

d = 2n'^ + n) 



Hp{k) = Fp{e^^^^''^V^We^' 



Kf ki 



+ V^YV) , (59) 



with p £ (IjCXj). Such a function is well defined for any 
channel and input covariance matrix. The form of the Ki is 
the one given in App. |X] By definition, M is the minimum 
of the function Hp{k) and, because of Eq. (I57> . Hp{0) — M. 
Furthermore, such a function is differentiable in fc = 0. There- 
fore, if the covariance matrix 7 is indeed optimal the function 
Hp{k) has to be stationary (critical) in fc = 0. This constraint 
is explicitly expressed by 



_d_ 

dki ^,^p 

for alH = 1 , . . 

_d_ 

dk, 



Fpie^U'^^^T^'e^U'^^K, 
d. We have that 

=0 
{Kfi + iK,) . 



y') - , (60) 



^'^Uk.K^ 



Y') 



(61) 



Since 7' + Y' is in Williamson form, we can apply the results 
of Apd.IbI On using Eq. ( I61> and Eq. (IB6> we obtain for the 
first derivative of the symplectic eigenvalues of the output 



'dk. 



tT,[K^-f' + ^'Ki_ 



(62) 



fc=0 



where tr^ denotes the trace of the leading principle submatrix 
corresponding to the mode j (see App.IbI. For the derivative 
of the output p-purity we hence obtain 



dki 



fe=0 



Fp(e'>^tik.KTy^EUk.K.^Y') 



(63) 



-Fp{j' + r) 






1 /(-) 



tiAKl^-f' + j'K, 



where /' = df /dx. In order for 7' to be the an input corre- 
sponding to a critical point, this derivative has to be zero for 
all generators Ki of Sp{2n, R). As can be promptly verified 
exploiting the form of the generators and Eq. ( I63> . the condi- 
tion of Eq. ( I60> results, for the symmetric Ki, in the following 
constraints on the submatrices of 7' ai, and Pki 



ails 



Pki^ 



h' 



kl 



h" 
"kl 



"kl "kl 



(64) 



for some real a, > 1 and 5'^,;, b'li G M. Since 7' is the covari- 
ance matrix of a pure state, all its symplectic eigenvalues are 
equal to 1, implying that \iay\ = l2n- Therefore 



-(77 (77 



l2n ■ 



Applying the previous condition to the submatrices of Eq. (|63 
yields 



l=i!^i 



U'? 



) = 1 



(66) 



which is equivalent to 



l,b',i = bl,=0. 



(67) 



Eq. J67t shows that the unique 7' that is consistent with 
a critical point is the identity l2n, corresponding to the 
?i-fold tensor product of a coherent state. It is easy to verify 
that the identity also satisfy Eq. ( I63> for the antisymmetric Ki. 



Summarizing, we have shown that the unique optimal input 
7 corresponds to 7' = V'^jV — 1, where 1/ is a symplectic 
transformation for which 



V^jV + V^YV 



(68) 



is in Williamson standard form. However, as the identity and 
Y themselves are in Williamson form, it is immediate to see 
that V = 1, yielding 7=1, which completes the proof. D 

For p = 2 the above multiplicativity result can easily be ex- 
tended to products of arbitrary channels acting on any number 
of modes without imposing additional constraints on the de- 
terminants of the Xi (apart from being non-zero): 

Proposition 3. - The maximal Gaussian output 2-purity of 
a tensor product of arbitrary multi-mode Gaussian channels 
$0Xi,eYi with det[Xi] ^ for all i is multiplicative. 

Proof. Because of Eq. j29> . the Gaussian multiplicativity issue 
reduces for p = 2 to the multiplicativity of the infimum of 
det[c/)x,y(7)] over all covariance matrices corresponding to 
pure Gaussian states. For a given Gaussian channel $x,f, 
making use of the Binet theorem (det AB = det A det B) 
and defining 

Y' = X-^YX-\ (69) 

one gets 

dct[X'^j{l,z)X + Y] :=^dct[Xfdet[^{l,z) + Y']. 

However, Y' can be diagonalized to Y' by a symplectic block 
matrix, which in turn does not change the determinant; 



Y' = SYS^ , ^-05, 



(70) 



Therefore, the problem is equivalent to verifying the multi- 
plicativity of the infimum of 



detm,z) + Y']^F2i^^y,)/r, 



(71) 



(65) which we know to hold true because of Proposition 2. D 



Note that Proposition 2 also implies multiplicativity for 
other multi-mode Gaussian channels. In particular if Xi = 
Xi Si are proportional to symplectic transformations with Xi > 
0, then multiplicativity holds for any $©Xi,eY; within the en- 
tire range p G (1,00). 



VI. COMMENTS AND OUTLOOK 



Acknowledgements 



We have addressed the multiplicativity of the maximal out- 
put p-purities of tensor products of Gaussian channels de- 
scribed by Eq. ( I18> . We have proved that, restricting to Gaus- 
sian inputs, the maximal output p-purities are multiplicative 
for any p G (l,oo) for single-mode channels with X = 
©"^iX,, Y = ®:[Li^i if dot Xt is the same for i = 1, ..., n, 
and that the ordinary 'purity', corresponding to p = 2, is 
multiplicative for any generic choice of multi-mode Gaussian 
channels. In particular, the maximal output purity is multi- 
plicative for identical n-fold single-mode channels for all val- 
ues of p € (1, oo) and in this case the optimal product output 
(which is independent of p) majorizes any other Gaussian out- 
put state. 

The restriction to Gaussian states, formally expressed by 
the definition of Gaussian maximal output purity of Eq. OH . 
is here motivated by essentially three arguments. Firstly, the 
question is interesting in its own right: Gaussian states have 
a prominent role in quantum information and communication 
with continuous variables, where many protocols completely 
rely on such states. Colloquially, one may say that the re- 
sults indicate that entangled Gaussian input states suffer more 
decoherence than uncorrected ones. The results and argu- 
ments presented in this paper constitute a strong hint towards 
the multiplicativity of the maximal output purities of general 
products of Gaussian channels under Gaussian inputs. In other 
words, it seems plausible that, in a fully Gaussian setting, in- 
put entanglement does not help to better preserve output pu- 
rity of quantum channels. Notice also that our proofs of mul- 
tiplicativity encompass several instances of interest. In par- 
ticular. Proposition 1 represents the case of the subsequent 
uses of a single channel, where input correlations could be 
distributed in time over the global input. In such an instance 
our result proves that Gaussian entangled input states suffer 
more decoherence than uncorrected states. We mention as 
well that Proposition 3 includes the relevant case of dissipa- 
tion of multi-mode systems in Gaussian reservoirs, provided 
that the coupling to the reservoir is the same for any mode, 
but allowing for generally different reservoir states in differ- 
ent modes [see Eq. ( I22H . 

Secondly, the maximal output purities under Gaussian in- 
puts readily deliver bounds for the maximal output purities 
of Gaussian channels, not restricting to Gaussian inputs. Fol- 
lowing the results presented in Refs. P9ll one might conjecture 
that these bounds are tight and that Gaussian input states are 
already optimal. 

Thirdly, the issues considered here could be the first steps 
towards a general theory of quantum information of Gaus- 
sian states, linking output purities to the Gaussian instance of 
the entanglement of formation |IT§], and Gaussian versions of 
channel capacities. As such, the Gaussian picture would de- 
liver a convenient and powerful testbed in entanglement the- 
ory, in an instance for which a complete solution for the seem- 
ingly unrelated additivity and multiplicativity problems may 
be anticipated. It is the aim for future work to establish this 
connection in generality. 
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APPENDIX A: GENERATORS OF THE SYMPLECTIC 
GROUP 

As can be easily verified from the expression of the condi- 
tion S'^crS = a, the symplectic group Sp{2n, R) is generated 
by those matrices which can be written as K = a J, where J 
is a symmetric 2n x 2n matrix I12ll . The antisymmetric gener- 
ators result in orthogonal symplectic transformations, giving 
rise to the compact subgroup K{n) = Sp{2n, R) n S0{2n). 
Such a subset of transformations is constituted by 'energy pre- 
serving' or passive operations. In contrast, the symmetric gen- 
erators generate the non compact subset of the group (made up 
of active transformations, like squeezings). A basis of such 
generators can be built by means of transformations affecting 
only 1 or 2 modes at a time. We define 



P-- 



1 
-1 



1 

1 



(Al) 



and recall the definition of a in Eq. (to be understood in 
the single mode, 2x2 instance). Single mode transformations 
are generated by 



-cr/3 = (5 , aS = P , 



(A2) 



where a generates the compact single mode rotations while /3 
and S generate single mode squeezings. 
Two-mode transformations (corresponding to the compact 
set) are generated by 



a 
a 



-1 

1 



(A3) 



Whereas two-mode transformations (corresponding to the 
non-compact set) are generated by 



5 
5 



p 
(3 



The complete set of generators 

{Ki : i = l,...,2n^ + n} 



(A4) 



(A5) 



is described by Eq. JA2l i for any mode and by Eqs. jA3t and 
( IA4> for any couple of modes. The total number of indepen- 
dent generators is 



3n + 4n(n - l)/2 = 2n'^ + n. 



(A6) 



The number of generators of the compact subgroup, which we 
refer to as {Li : i == 1, ..., n^} in this paper, is n + 2n{n — 
l)/2 = n\ 



APPENDIX B: SYMPLECTIC PERTURBATIONS 



with eigenvectors 



We consider a covariance matrix of n modes in Williamson 
form 



(Bl) 



and investigate the variations of the symplectic eigenvalues 
vi, ...,v^ under an additive perturbation. Let us consider 7 + 
kP, with k £ ]R+, and P being a symmetric 2n x 2n matrix, 
partitioned in terms of 2 x 2 submatrices Pjj as 




/ Pii P12 



Pin \ 



P'^ 






P^ 



^n — 1 n 
P 



The eigenvalues of the matrix 10-7 are given by 



(B2) 



(B3) 



Vj+ = (0,...,0, z,l, 0,...,0)^ (B4) 

mode j 

V,. = (0,...,0, l,i, 0,...,0f, (B5) 

mode j 

for j = 1, . . . , n, SO that ia^Vj^ = ^u^Vj^. Now, one has 



d 
dk 



Vj+{i(TP)vJ_^ ^Vj + P' 



fc=0 



"j+ 



trP,-, 



tr,P , 



(B6) 



where we have defined tr^ as the trace of the leading subma- 
trix associated to mode j. The first order derivative of the 
symplectic eigenvalue v^ is just given by the trace of the 2x2 
principal submatrix related to mode j of the matrix embody- 
ing the perturbation. 
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